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Abstract. - We perform molecular dynamics simulations of a model glass-forming liquid to 
measure the size of kinetic heterogeneities, using a dynamic susceptibility Xss(a,t) that quantifies 
the number of particles whose dynamics are correlated on the length scale a and time scale t. By 
measuring Xmiflji) as a function of both a and t, we locate local maxima \* a t distances a* and 
times t*. Near the dynamical glass transition, we find two types of maxima, both correlated with 
crossovers in the dynamical behavior: a smaller maximum corresponding to the crossover from 
ballistic to sub-diffusive motion, and a larger maximum corresponding to the crossover from sub- 
diffusive to diffusive motion. Our results indicate that kinetic heterogeneities are not necessarily 
signatures of an impending glass or jamming transition. 



Introduction. — The dynamics in a liquid become in- 
creasingly correlated in space and heterogeneous in time 
as the system is cooled toward its glass transition. These 
correlations are most pronounced at time scales on the or- 

' der of the relaxation time of the system and are typically 
quantified using a dynamic susceptibility that measures 
fluctuations in the number of mobile particles [IH2]- A 
particle is deemed mobile if it has moved a distance of at 
least a in a time interval t, so the dynamic susceptibility 
depends on the choice of a and t. For a fixed at a fraction 

' of the particle diameter, the dynamic susceptibility ex- 
hibits a maximum at times t on the order of the relaxation 
time. This maximum increases in height as temperature 
is lowered towards the glass transition and the relaxation 
time increases [6ll8Hl8j. Thus, the peak in the dynamic 
susceptibility is viewed as an important signature of the 
glass transition. 

In molecular dynamics simulations of glass-forming liq- 
uids, the dynamic susceptibility is typically measured as a 
function of t with the length scale a fixed. Studies of the 
length-scale dependence have been more rare [T5I[ToHT5I 

119) . However, recent experiments on granular systems [20 
suggest that it is instructive to study the dynamic suscep- 
tibility as a function of a as well as t. This approach has 
the advantage of avoiding any arbitrariness in the choice 
of a or t. 

In this paper, we use the dynamic susceptibility to mea- 
sure the spatial extent of kinetic heterogeneities as a func- 
tion of distance a and time f in a model glass-forming 



liquid. At temperatures where there is a well-defined sub- 
diffusive plateau in the root mean square displacement, 
we find that there are iwo distinct local maxima in the 
dynamic susceptibility Xss(a, t) quantifying the spatial ex- 
tent of kinetic heterogeneities. These two maxima corre- 
spond to crossovers in the dynamics. The larger of the two 
maxima, which has been observed in previous simulations, 
occurs at late times at the crossover from sub-diffusive to 
diffusive motion. This maximum quantifies the strong and 
well-studied correlations in the dynamics that arise at the 
scale of the relaxation time. In addition, we find a sec- 
ondary maximum at much earlier times, corresponding to 
the crossover from ballistic to sub-diffusive motion. The 
presence of this second maximum indicates that the dy- 
namics also exhibit significant spatial correlations at time 
scales corresponding to the trapping of a particle in the 
cage of its neighbors. The discovery of a second maximum 
implies that kinetic heterogeneities are not just associated 
with the onset of slow dynamics near a glass or jamming 
transition, and that they should be regarded more gener- 
ally as symptoms of crossovers in the dynamics. 

Model. — Our model consists of disks interacting in 
two dimensions via the purely repulsive pairwise potential 

for ry > (Ti + (Tj , 

where e is the characteristic energy scale and Oi is the 
radius of the ith disk. We use equal mixtures of disks 
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of radii a and lAa and equal mass m. We choose units 
such that 2a, e, m, and ks set equal to 1. This sets the 
time unit 2aymje = 1. To solve the dynamics, we per- 
form molecular dynamics simulations at fixed temperature 
and pressure in square periodic cells containing N = 100, 
400, or 1600 disks. We employ Gaussian constraints [3T] 
to fix the instantaneous temperature T = \pi\ 2 /2(N — 1) 
and hydrostatic pressure p = (fi ■ Fi/2 + NT)/L 2 , where 
Fi = —VJ2j ^( r «i)j allowing the side length L of the pe- 
riodic cell to vary. In our runs, we fix p = 10~ 2 and 
measure the dynamics at various temperatures T = 10 -4 
to T = 9 x 10~ 4 by running 10 to 20 simulations for dura- 
tions between At = 10 5 and At = 10 7 after equilibrating 
at the temperature for one fifth as long. We arrive at the 
temperature by quenching at rates 10~ 9 or 10 -8 starting 
from well-equilibrated configurations at T = 10~ 3 . Re- 
sults here are shown for a quench rate of 10~ 9 . 

Overlap function and its correlations. Correla- 
tions in the dynamics are measured by using the standard 
time- and distance-dependent order parameter 

q s (a,t;i,r) = w a (fi(T) -^(r + t)), (2) 

where w a (f) is an overlap function |12) . Here, w a (f) = 1 
if |r| < a and otherwise. We denote the average over 
all particles i by Q s (a,t;T) = (1/^V) q s (a, t; i, r) and 
its time average by Q s (a, t) = (Q s (a, t; t)) t . At very short 
times, Q s (a, t) = 1 because no particles have moved a dis- 
tance a. At very long times, positions at time r + t are mi- 
correlated with positions at time r, so lim^oo Q s (a,t) ~ 
<fia 2 /N, where </> is the packing fraction. Between these 
two extremes, Q s (a,t) decays smoothly from 1 to (f>a 2 /N. 

The spatial extent of fluctuations can be characterized 
by the variance of the overlap function over different start- 
ing times [T2] 

Xss (a, t) = N ((Q s (a, t; t) 2 ) t - Q s (a, t) 2 ) . (3) 

As shown by the counting argument of ref. [19] . x ss (a,i) 
is roughly the number of particles that move a distance a 
over a time t in a correlated manner. 

Note that the spatially averaged overlap function 
Q s (a,t;r) is the self part of the order parameter 
Q(a,t;r) = (1/N) ^ io (fi(r) -fj (r+t)), and Xss(a,t) is 
the self-self part of the dynamic susceptibility Xi{ a ^) = 
N((Q(a,t;T) 2 ) T - {(Q(a, t; r)) r ) 2 ). Glotzer et al showed 
that the self-self part of the dynamic susceptibility domi- 
nates over the distinct-distinct and the self-distinct parts 
for a model glass-forming liquid [12II13] . 

Results. — The color plots of fig. [T] display the spa- 
tial and temporal dependence of the dynamic susceptibil- 
ity Xss(a, t)- The three plots represent three temperatures: 
one well above the dynamic glass transition T g 5 x 10~ 4 , 
one slightly above it, and one below T g . In each plot, 
we also display the root-mean-square (rms) displacement, 
Ar(t) = -J (\ri(0) — fi(t)\ 2 ). Not surprisingly, for a given t 
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Fig. 1: Color plots of the value of the dynamic susceptibility 
Xss(a, t) as a function of the lag time, t and the overlap dis- 
tance, a. Results are shown at a fixed pressure, p = 1(T 2 , and 
system size, iV = 1600, at three different temperatures: (a) 
T = 7 x 10 -4 , (b) T = 5.5 x 10~ 4 , and (c) T = 2 x 10~ 4 . The 
dynamical glass transition is at T ~ 5 x 10~ 4 . For compar- 
ison, we plot the root-mean-square displacement Ar(t) (solid 
curves). The circles mark the locations, (a*,C), of observed 
local maxima of Xas(a, i). Filled circles represent cage-escaping 
local maxima near the end of the sub-diffusive plateau of Ar(t). 
Such maxima are observed for T > 5 x 10 -4 . Open circles rep- 
resent cage-exploring local maxima near the beginning of the 
plateau. Such maxima are observed for T < 5.5 x 10~ 4 , so that 
there is a range of temperatures over which both local maxima 
exist. 



we find that Xss{a,t) is largest for a near Ar(t), the char- 
acteristic distance over which some, but not all, particles 
have moved. This is consistent with granular experiments 
by Lechenault, et al [2D] . We place circles on the locations 
(i*, a*) of the local maxima of Xss(a, t). 

For the highest temperature (fig. [Ha)), well above the 
glass transition, the rms displacement Ar(t) shows a bal- 
listic regime at small i, the hint of a sub-diffusive regime at 
intermediate times, and a diffusive regime at long times. 
We observe one local maximum of Xss(o ? i), marked by 
a solid circle, at the onset of the diffusive regime: the 
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heterogeneities are largest when some particles have be- 
gun to diffuse but others remain caged. Note also that 
the maximum at (t*,a*) lies somewhat below the Ar(t) 
curve, indicating that larger clusters of dynamically cor- 
related particles are formed by less mobile particles Q|. 

For the intermediate temperature (fig. IHb)), the sub- 
diffusive plateau of Ar(i) is well established and spreads 
out over two decades. Again, we observe a local maximum 
of Xss(o, t) near the onset of the diffusive regime. However, 
in addition to this "cage-escaping" maximum we also ob- 
serve a secondary "cage-exploring" maximum, marked by 
an open circle, near the onset of the plateau; the dynam- 
ics are also heterogeneous on time scales when particles 
are becoming constrained by their neighbors. Thus, in 
the regime where there are two distinct crossovers in the 
dynamics-one from ballistic to sub-diffusive motion, and 
one from sub-diffusive to diffusive motion-there are also 
two maxima in x ss (a,i), located at distances a and times 
t corresponding to the two crossovers. 

For the lowest temperature shown (fig. [2c)), we do not 
reach a diffusive regime on the time scale of our simula- 
tion, which means that this temperature lies below the 
dynamic glass transition temperature T g w 5 x 10~ 4 . We 
therefore do not observe the cage-escaping maximum as- 
sociated with the onset of the diffusive regime. However, 
the location of the cage-exploring maximum remains rela- 
tively unchanged, as shown by the open circle. Although 
the system is out of equilibrium for all T < T gi we observe 
mild aging effects only for the highest of these tempera- 
tures studied, T = 4 x 10 -4 and T = 4.5 x 1CT 4 , where the 
equilibration times are comparable to the quench times, 
waiting times, and run times of our simulations. 

To study the evolution of the maxima with tempera- 
ture, we plot in Fig. [5] Xss(i, as a function of t for many 
different fixed values of a covering the same range of a 
as in Fig. [1] The thick red curve marks the maximum 
of Xss with respect to a at each t. Maxima of this curve 
(marked with circles) therefore represent maxima of x ss 
with respect to both a and t. Fig. [2] shows that at high 
temperatures, only one local maximum of x ss can be re- 
solved. This is the well-known maximum at times t com- 
parable to the relaxation time; it is represented by a filled 
circle in fig. [21(a). As T is lowered, however, a shoulder 
becomes apparent at time scales corresponding to cage 
exploring, as shown in fig. (2fb). This shoulder evolves 
into a secondary maximum, denoted by an open circle in 
fig. [U(c), as T is lowered still further. Note that the pri- 
mary cage-escaping maximum is quite broad in t and that 
the secondary cage-exploring maximum is only observed 
at temperatures at which the cage-escaping maximum has 
moved to sufficiently long times. Thus, the reason why the 



1 In all cases, we also find that Xm (i, t) increases along Ar(t) 
with decreasing t for the shortest times observed, well within the 
ballistic regime, reflecting correlations in the velocities. We observe 
this behavior at all temperatures. This is an artifact of the barostat 
used; we do not observe this increase when the system is held at 
fixed area. 
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Fig. 2: Dynamic susceptibility Xss(a,t) vs lag time t for p = 
10 -2 , N — 1600, and four different temperatures above but 
near the glass transition: (a) T = 6.5 x 10~ 4 , (b) T = 6 x 
10~ 4 , (c) T = 5.8 x 10" 4 , and (d) T = 5.5 x 10~ 4 . At each 
temperature, we show multiple curves, each corresponding to a 
different value of the overlap distance a. The range of overlap 
distances represented by the curves is the same as in the color 
plots of fig. [T] The maximum of Xss with respect to a at each 
value of t is represented by the heavy red curve. Circles mark 
the lag time t* and amplitude x* °f local maxima of Xss(a,t). 
As in fig. [U filled circles represent cage-breaking local maxima 
near the end of the sub-diffusive plateau of Ar(t), and open 
circles represent cage-forming local maxima near the beginning 
of the plateau. 



observation of two distinct maxima is restricted to tem- 
peratures fairly close to the glass transition is because the 
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Fig. 3: Local maxima of the dynamic susceptibility. We plot 
(a) the magnitudes of all observed local maxima of Xss(o, t), 
X* , (b) their lag times t*, and (c) their overlap distances a* 
as a function of temperature at fixed pressure p = 10 -2 and 
system size JV = 1600. As in fig. [T] the filled circles represent 
cage-escaping maxima, while the open circles represent cage- 
exploring maxima. The vertical dashed lines denote the range 
over which we observe both types of maxima. Note that (a) and 
(c) are presented on a linear scale, while (b) is presented on a 
log-linear scale. Error bars represent the larger of (1) the stan- 
dard deviation of the mean calculated from simulations with 
different initial conditions and (2) the uncertainty associated 
with the discrete sampling of a and t. 



smaller cage-exploring maximum is obscured by the larger 
cage-escaping maximum at higher temperatures. 

Fig. [3] shows the temperature dependence of the two 
maxima of the dynamic susceptibility. We plot the ampli- 
tude of the maxima x* an d their associated locations in 
time and space, t*, and a* , as functions of T at fixed pres- 
sure p = 10~ 2 and system size N — 1600. The amplitude 
X* and lag time t* of the primary, cage-escaping maxi- 
mum (solid symbols) both increase as T decreases until 
t* passes beyond our observable time window, in accord 
with many previous observations. The overlap distance 
a* of the cage-escaping maximum is on the order of half a 
disk diameter for temperatures near the glass transition, 
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Fig. 4: Dependence of the relaxation time t* on the spatial 
extent of kinetic heterogeneities, x* ■ We plot the position of 
the cage-escaping maximum, t* vs the magnitude of that maxi- 
mum, x* a t fixed pressure p = 10~ 2 . Each point corresponds to 
a different temperature above T g . Results are shown for three 
different system sizes N, as labelled, on both (a) a log-linear 
scale and (b) a log-log scale. The straight lines in (a) and (b) 
are exponential and power law fits, respectively. Error bars 
represent the larger of (1) the standard deviation of the mean 
calculated from simulations with different initial conditions and 
(2) the uncertainty associated with the discrete sampling of a 
and t. 



a distance about two times larger than the plateau of the 
rms displacement. This distance presumably reflects the 
scale of the rearrangements necessary for diffusive motion 
to occur. At high temperatures above T = 6 x 10 -4 , the 
plateau in the rms displacement is less well defined and 
the cage-escaping maximum lies at slightly larger values 
of a*. 

The temperature dependence of the cage-exploring, sec- 
ondary maximum is markedly different from that of the 
primary maximum. Figure [3] shows that x* an d t* re- 
main roughly constant as temperature decreases from 
T = 5.5 x 10~ 4 , just above T g w 5 x 10-4, to T = 10" 4 , 
which is well below it. This is consistent with the behav- 
ior of the rms displacement; the time of the onset of the 
plateau does not vary much with temperature as long as 
the temperature is low enough for the plateau to exist. 
The cage-exploring overlap distance a* decreases with de- 
creasing T, remaining near the plateau value in the rms 
displacement. 
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The amplitude x* of the cage-escaping maximum is a 
measure of the number of disks whose dynamics are cor- 
related at the crossover to diffusion, while t* is a measure 
of the timescale for this longest, or a-, relaxation process. 
In fig. 2] we plot x* i the number of disks whose dynam- 
ics are correlated at the cage-escaping maximum of x ss , 
as a function of the time scale corresponding to the cage- 
escaping maximum, t* . This time scale is a measure of the 
a-relaxation process. We have plotted x* vs. t* for three 
different system sizes, N — 100, 400, and 1600. Note 
the substantial system-size dependence in x* [22], even 
for temperatures for which the size of correlated regions 
n* rs Ax* is only on the order of 10. This indicates that the 
finite size of the system affects the dynamics even when the 
relative linear size of the correlated regions is as small as 
y/n*/N « 0.1. This strong effect is likely due to a broad 
distribution of sizes of correlated regions. We find that 
while the spatial correlation function of the overlap func- 
tion has an average correlation length consistent with the 
value \J ir(<r 2 )n* /<fi expected from a compact correlated re- 
gion, it also has a long-ranged tail, decaying slower than 
exponentially with distance. This long range tail reflects 
intermittent periods when the the clusters are extended. 
We find that t* and a* also depend on system size. If a 
is held fixed, t* decreases and x* increases with increas- 
ing system size N at all temperatures studied, consistent 
with Ref. {22] . However, when we maximize over a, we find 
that at sufficiently high T, t* increases with N because a* 
increases with N there. 

Fig. 2] shows that for each system size N, the depen- 
dence of x* on t* is qualitatively similar. Note that the 
log-linear plot (fig. Ufa)) is straighter than the log- log plot 
(fig- HP 3 )) for each system size, indicating that the depen- 
dence of t* on x* is described better by an exponential 
than by a power law. Such an exponential dependence of 
relaxation time on cluster size is expected in scenarios in- 
volving cooperative rearrangements |14) . Previous numer- 
ical studies of model glass formers fit x* an d t* to power 
laws in T—T c , where T c is a critical temperature [12|13|15] . 
implying that the relationship between x* an d t* is also 
a power law. However, these studies determined x* as 
the maximum of Xi{ a i^) with respect to t at fixed a, in- 
stead of searching for a maximum with respect to both t 
and a. Also, the relationship between x* an d t* was not 
shown directly. An analysis of a lower bound of Xi f° r 
several real glass-forming liquids found an exponential de- 
pendence near the glass transition [5] , consistent with our 
results. 

Discussion. — We quantified the size of kinetic het- 
erogeneities for a model glass-forming liquid as a function 
of temperature in terms of fluctuations of a dynamic or- 
der parameter. By measuring the dynamic susceptibility 
Xss(a,t) as a function of both lag time and overlap dis- 
tance, we identified two distinct local maxima. The pri- 
mary, "cage-escaping" maximum, which has been found 
in previous numerical and experimental studies, occurs at 



the end of the plateau in the root-mean-square displace- 
ment, near the onset of diffusive transport. We also dis- 
covered a secondary, "cage-exploring" maximum at the 
beginning of the plateau in the mean-squared displace- 
ment, which corresponds to the ballistic to sub-diffusive 
crossover. Thus, the dynamics are most heterogeneous 
at displacements and times corresponding to crossovers in 
the dynamics. 

The secondary, cage-exploring maximum at the 
crossover from ballistic to sub-diffusive motion emerges 
from the shoulder of the primary cage-escaping maximum 
when the temperature is low enough that there is a well- 
defined plateau in the mean-squared displacement. We 
find that the plateau must extend over at least two orders 
of magnitude in order for the cage-exploring maximum to 
appear. We observe both maxima over the range of tem- 
peratures that are low enough that the mean-squared dis- 
placement exhibits a broad plateau but high enough that 
the eventual crossover to the diffusive regime is still ob- 
servable. Note that although both maxima emerge, grow, 
and shift continuously with temperature, there is an ap- 
parent discontinuous drop in the largest observed maxi- 
mum at the dynamical glass transition. This occurs when 
the primary maximum passes out of observation range, 
leaving only the secondary, cage-exploring maximum (see 

fig. Bp))). 

While the crossover from sub-diffusive to diffusive mo- 
tion measures the approach to the glass transition, the 
crossover from ballistic to sub-diffusive motion is a general 
and relatively innocuous feature of dense fluids and solids. 
Our finding that both crossovers exhibit maxima in the dy- 
namic susceptibility suggests that spatially heterogeneous 
dynamics are a general feature of dynamical crossovers 
caused by interactions. Indeed, a crystalline system of 
monodisperse disks exhibits a cage-exploring maximum 
at the end of the ballistic regime that is similar to the 
one shown here. We suspect that more exotic dynamical 
crossovers, like those between two diffusive regimes with 
different diffusion coefficients, should also exhibit kinetic 
heterogeneities. The existence of kinetic heterogeneities 
is therefore a general feature of dynamical crossovers and 
should not necessarily be taken as a sign of an impending 
glass transition. 
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